Visualization of vortex bound states in polarized Fermi gases at unitarity 
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We theoretically analyze a single vortex in a spin polarized 3D trapped atomic Fermi gas near a 
broad Feshbach resonance. Above a critical polarization the Andreev-like bound states inside the 
core become occupied by the majority spin component. As a result, the local density difference 
at the core center suddenly rises at low temperatures. This provides a way to visualize the lowest 
bound state using phase-contrast imaging. As the polarization increases, the core expands gradually 
and the energy of the lowest bound state decreases. 
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The achievement of superfluidity in trapped ultra-cold 
atomic 6 Li gases is a landmark advance in the history of 
physics [H- This is attained by utilizing a broad Fesh- 
bach resonance, which is used to tune the inter-atomic 
interactions. By changing the inverse scattering length 
a s continuously from negative to positive values, a two- 
component Fermi gas with equal spin populations has 
a ground state which crosses smoothly from Bardeen- 
Cooper-Schrieffer (BCS) superfluidity to a Bose-Einstein 
condensate (BEC) of tightly bound pairs. Of particu- 
lar interest is the unitarity regime near resonance, where 
the scattering length diverges (l/a s — 0). Since the 
inter-particle spacing is the only relevant length scale, 
the Fermi gas exhibits a universal behavior 0|. 

Quantized vortices are a clear-cut confirmation of su- 
perfluidity, and were demonstrated experimentally by 
Zwierlein et at The equilibrium properties of 

vortices in a symmetric Fermi superfluid at crossover 
have been the subject of intense theoretical studies 
@, 0, @, IE 0, S, H- The Andreev-like bound states, 
which are the fermionic quasiparticle excitations local- 
ized in the core, have been widely discussed (!, IEBI11- 
These bound states are found to play a key role in the 
structure of vortices. 

Most recently, Fermi gases with unequal spin popula- 
tions have been the subject of considerable experimental 
13, [LH and theoretical interest [H, El Ei EE EE E3, El • 
The presence of spin polarization leads to exotic forms of 
pairing, such as breached pairing [12j or Sarma super- 
fluidity [ijj], phase separation [lj|, and spatially mod- 
ulated Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states 
An agreement on the true ground state of polar- 
ized fermionic superfluidity is yet to be reached. How- 
ever, three recent measurements on the density profiles of 
polarized 6 Li gases [IE EH, near a Feshbach resonance, 
indicate a paired superfluid core surrounded by the ex- 
cess unpaired fermions consistent with a picture of phase 
separation. 

Combining spin polarization with a vortex may help to 
resolve the issue of the nature of polarized fermion pair- 



ing. It is natural to ask how unequal spin populations 
affect the vortex structure, and how vortex bound states 
evolve as the polarization increases. This issue arises in 
the context of pairing and superfluidity in many fields of 
physics [ijj. It is highly relevant to the condensed mat- 
ter community, where polarized superfluidity is created 
by applying a magnetic field. There is now strong exper- 
imental evidence for the existence of FFLO states in the 
heavy fermion superconductor CeCoIns under high fields 
[icj ]. Strongly interacting polarized Fermi gases have also 
been under close scrutiny in nuclear matter |2l fl. neutron 
and high density quark matter 



stars 
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where 



the spin polarization is created by differences between 
chemical potentials and/or by mass differences between 
fermions that form pairs. Polarized vortices of color su- 
perfluidity in rotating neutron stars are a possible mech- 
anism for observed glitches in pulsar timing [lit ]. 

Here we investigate the properties of a singly quantized 
vortex in polarized atomic Fermi gases at unitarity, in a 
cylindrically symmetric trap. Our main results are: 

(A) We clarify the density profiles of both spin compo- 
nents as a function of polarization. In addition to phase 
separation, the vortex core suddenly accommodates the 
excess majority fermions above a critical polarization or a 
critical chemical potential difference, resulting in a rapid 
rise of the local density difference inside the core. 

(B) The local fermionic density of states explains the 
sudden appearance of an unpaired core of excess major- 
ity atoms at the vortex center. The Andreev-like bound 
states in the core are occupied when the critical chemi- 
cal potential difference equals the lowest available energy. 
This provides a clear visualization of vortex bound states 
using phase-contrast imaging [24]. 

(C) With increasing polarization, the vortex core ex- 
pands while the lowest bound state energy decreases. 

The above results are obtained by numerically solving 
the mean-field Bogoliubov-de Gennes (BdG) equations in 
a fully self-consistent fashion [E E3], assuming a pairing 
order parameter that preserves the cylindrical and axi- 
ally translational symmetries. Symmetry breaking is also 
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possible, i.e., the order parameter may distort cylindri- 
cally. This scenario merits further study. 

Fermi gases of 6 Li atoms near a broad Feshbach reso- 
nance are well characterized using a single channel model 
[25l | . The BdG equations describing the quasiparticle 
wave functions u v (r) and v v (r) , with excitation energies 
E v read 0: 



Ho - m 
A*(r) 



where Ho 
muj 2 (x 2 



A(r) 
-H + hi 



v n (r) 



E„ 



u n (r) 
Vr, (r) 



(1) 



-h 2 V 2 /2m + V ext (r), and V ext (r) 



yj/2 is the transverse trapping potential. 
Along the z axis we instead assume free motion over a 
length L. To account for the unequal spin population 
N a for a =t)l) the chemical potentials are shifted as 
/i^l = p ± leading to different quasiparticle wave 
functions for the two components. However, there is a 
symmetry of the BdG equations under the replacement 

u nl ( r ) ~* v rt ( r )' v vi ( r ) ~* ~ U ')T ( r ); E vl -> ~ E vt- We 
can thus retain only (r) and u^f (r) in Eq. (JTJ) , and 

keep solutions with both positive and negative energies. 

The order parameter A(r) and the chemical potentials 

/if 4 are determined by self-consistency equations for the 

gap, A(r) = gJ2 v u v( r ) v v( r ) f ( E v)i and the particle 

density of each component: nf (r) = J2 V \ u n ( r )| 2 / ( E v) 

and n x (r) = l«n (r)| 2 / (-£„) 

constrained so that J" rfrrt CT (r) = iV, 



,x/k B T 



These must be 
„ , where / (x) 



■ l) is the Fermi distribution function, and 



g (< 0) is the bare coupling constant, which is related to 
the s-wave scattering length via the regularization pre- 
scription: (47rS 2 o s /m) = 1/g + J2k l/2ek- 

We solve these equations via a hybrid procedure, 
by introducing a high energy cut-off E c above which 
we use a local density approximation (LDA) for high- 
lying excitation levels. The standard regularization 
prescription then yields an effective coupling con- 
stant through the self-consistency equation A(r) = 
9eff ( r ) Y^n Ur i ( r ) v n ( r ) / i E v)> wnere the cut-off summa- 
tion Y^' v is now restricted to \E V \ < E c . Further details of 
this will be given elsewhere. A clear limitation of the pro- 
cedure is the use of mean-field factorizations implicit in 
the BdG equations. From earlier work, we expect this to 
neglect quantum fluctuations that alter the ground-state 
energy, while remaining qualitatively correct [27j |. 

Below the cut-off, we solve the BdG equations by work- 
ing in cylindrical coordinates (p, if, z) and taking A(r) = 
A(p)e~'" p for a singly quantized vortex. Assuming peri- 
odic boundary conditions at z — ±L/2, we write, for the 
normalized modes, u v (r) = u nm k^ (p) e lmv e lk * z / \j2nL 
and v n (r) = v nmkss (p) e^ m+1 ^e ik - z /V2nL with k z = 
2irl/L. As a consequence, the BdG equations decouple 
into different m and I sectors 0. Expanding the radial 
functions u„ m k z (p) and v nm k z (p) in a basis set of 2D 
harmonic oscillators, we then solve a matrix eigenvalue 
problem in each sector. 




Figure 1: (Color online). Density profiles of the ma- 
jority (j-state, solid lines) and minority (J.-state, dashed 
lines) components at T = 0.05IV for three typical value 
of polarization: p = 0.12 (a), p = 0.35 (b), and 
p = 0.75 (c). Density differences are also plotted in 
dotted-dashed lines. All the profiles are normalized by 
n a ,TF= (1 + f3y 3/s v / 15tvNX/2/ (6tt 2 ) {h/mtu)- 3 ' 2 , which 
is the peak density for a symmetric gas at unitarity. Panel 
(d) shows the order parameter profiles. The small oscillations 
at the edge are a finite size effect. 



In greater detail, we consider a gas at unitarity with 
the number of total atoms in the range N = + N± — 
2 x 10 3 — 4 x 10 4 . Two characteristic scales may be defined 
by considering a symmetric ideal Fermi gas at zero tem- 
perature. In the LDA analysis this leads to a Thomas- 
Fermi (TF) radius p° TF = (157rNA/2) 1/6 ^h/muj, and a 

1/3 

Fermi energy E F = (157rAfA/16) i/J nw = k B T F , where 
we define A = Lj p% F as the aspect ratio of the trap. 
Throughout this Letter, we calculate results at the Fes- 
hbach resonance with l/a s = and use A = 1 and 
E c ~ 2Ef- We also considered coupling constants in 
the BCS regime but observed no significant changes. Di- 
mensionality effects will be treated elsewhere. 

Numerical accuracy was checked by increasing E c up 
to AEp. Due to the high accuracy of our hybrid cut-off 
procedure, the results were found to be essentially in- 
dependent of the cut-off energy. We note finally that, 
for a symmetric gas at unitarity, universality implies a 
TF radius of ptf — (1 + p% F , a chemical poten- 

N3/5 



tial p = (1 + 0) Ef, and a maximum order parameter 
A = 8 (1 + (3) 3/5 E F /e 2 0], where BCS theory predicts 
the universal parameter (3 ~ —0.41. 

We present in Figs, la, lb and lc the density profile of 
each component, as well as the density difference Sn (r) = 
7if (r)— (r), for several polarizations p = (JVf — N±) /N 
at T = 0.05Tp and TV = 10 4 . Because of the uniform dis- 
tribution along the z axis, these profiles are linked to 
the experimentally observed column densities in the ax- 
ial direction. Apart from the apparent phase separation 
at the edge, the most salient feature of the figures is the 



Figure 2: (Color online). Left panel: centre density differ- 
ence as a function of polarization at N = 10 4 . Inset shows 
the dependence on the chemical potential difference. Right 
panel: critical polarization and critical chemical difference as 
a function of the number of total particles. 



Figure 3: (Color online). Local fermionic density of states of 
spin up (solid lines) and spin down (dashed lines) components 
inside the vortex core at N = 10 4 and T = 0.05Tf. The thin 
line in (a) shows the LDOS at p = 0. Arrows points to the 
position at the effective energy of the lowest bound state. 



development of a polarized normal shell inside the vortex 
core above a certain polarization. This is clearly visible 
as a prominent peak in the density difference, of width 
about 0.05ptf • This is observable in the column inte- 
grated density difference, which is directly measurable 
by phase-contrast imaging lfj, 24 1. 

The onset of a polarized normal shell at the core cen- 
ter is demonstrated by the central density difference as 
a function of the polarization. This is shown in Fig. 2a, 
which represents the most important result of this Let- 
ter. At a sufficiently low temperature, i.e., T = O.OlTp, 
a sudden rise of the center density difference appears at 
a critical polarization p c ~ 0.30. The critical chemical 
potential difference is Sp c ~ 0.36-Ef ~ Aq/2Ejp, with a 
transition width of around ksT. This transition is there- 
fore much smoother at finite temperature. The critical 
polarization is nearly independent of the overall number 
of atoms N, as shown in Fig. 2b for N up to 4 x 10 4 . 
We therefore expect that this will apply to current ex- 
periments, where the typical number of atoms is around 
10 5 , and would survive even in the thermodynamic limit. 

The appearance of this intriguing shell structure is 
closely related to the Andreev-like bound states inside 
the core. In the BCS regime, these states are formed by 
the spatial variation of the order parameter around the 
center (see, i.e., Fig. Id), analogous to a potential well 
for quasiparticles, of depth Ao and of radius equal to the 
coherence length £ = Hvf /Ao. Hence, the confinement of 
the well gives rise to discrete bound levels with spacing 
of order h 2 /mt; 2 = A\j2E F 0|. This qualitative picture 
persists in the strongly interacting unitarity limit [§]. 

To provide an intuitive explanation of our results, we 
calculate the local density of states (LDOS), 

Ni(r,E) = J2 \u n (v)\ 2 5{E-E n ), 
N^E) = J2 K(r)\ 2 HE + E v ) . (2) 

At low temperature, when integrated over negative en- 
ergy, this leads to the density profiles n a (r). In Fig. 3 



we show how the LDOS inside the core evolves with in- 
creasing the polarization. A small spectral broadening 
of about 0.01-Ef has been used to regularize the delta 
function. Without any polarization the LDOS of the two 
components coincides, leading to a sharp peak located 
at positive energy E® s ~ Ag/2Ep, associated with the 
lowest Andreev-like bound state. 

In the presence of spin-polarization the peak in the 
density of states shifts in different directions for the two 
components. To a good approximation, the energy sepa- 
ration between the two peaks at the vortex center equals 
28(i. Thus, in the general case of a nonzero polarization 
one may define an effective energy of the lowest bound 
state, Eb s , as the midpoint of these two peaks located 
at Eb s =F 8(i. Therefore, a net density difference results 
precisely when the peak in JVj (r = 0, E) crosses zero en- 
ergy i.e., 8(i = Eb s . This results in a bound state for 
the majority spin component, which explains why a po- 
larized normal shell emerges above a critical population 
chemical potential 8fi c ~ E® s ~ A\j2Ep. 

Thus, the integrated column density difference is an 
indicator of the lowest vortex bound state, and a mea- 
surement of the critical polarization p c gives its energy. 

We now consider the dependence of the vortex core size 
on the polarization. We extract the core size from the 
superfluid density n s (r) , defined as a ratio of the current 
density j (rl = n s (r) v s to the superfluid velocity v s = 
(H/2mp) f> [9(|, where, since our normal fluid solutions are 
non-rotating: 

jW = ^E K d <P u nf ( E v) + "v^f (~ E n)] 0- (3) 
mp ^ 

The resulting superfluid density profiles are plotted in 
the inset of Fig. 4a. The core size may be quantified as 
the distance from the vortex core at which the superfluid 
density is 90% of its maximum value, namely, £90. From 
Fig. 4a, the core size increases gradually with increasing 
polarization, and almost doubles at large polarization. 

To explain this, note that while a phase separation 
occurs at any nonzero polarization, only the unpolar- 
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Figure 4: (Color online). Vortex core size (a) and the effective 
energy of the lowest bound state (b) as a function of polariza- 
tion at N = 10 4 and T = 0.05Tf. The core size £go at p = 
is about 2.5k p 1 , where Jcf non-interacting Fermi wavelength 
at center. Solid lines are the scaling relations as described in 
the text. Inset shows superfluidity density profiles. 

ized superfluid part can form a vortex. Thus, the vortex 
core should expand with a scaling of £90 oc (2iVj) cx 
(1 — p) In Fig 4a this scaling is plotted by a 

solid line, which fits well with our numerical results. Ac- 
cordingly, one may suspect that the energy of the lowest 
bound state will decrease as Et, s oc 1/ £f oc (1 — p) 2 ^ 3 . 
This is consistent with the effective energy of the lowest 
bound state shown in Fig. 4b. We expect a phase sepa- 
ration into multiple vortex cores in a vortex lattice, as in 
current non-polarized experiments []]]. 

We have considered an aspect ratio A = 1, which is 
closer to the MIT experimental setup [l(J than the Rice 
experiment (which has A = 50 In the opposite limit 

of A <C 1, an interesting aspect of dimensionality would 
arise. Due to strong phase fluctuations, this quasi-2D 
geometry would favor the spontaneous formation of vor- 
tices at finite temperature [26]. As a result, a lattice 
of vortex-anti- vortex pairs without phase separation may 
emerge as the ground state. In such a configuration, the 
spin polarization would be sustained by a polarized nor- 
mal shell inside the vortex cores, analogous to a type-II 
superconductor in a magnetic field. 

In conclusion, we have analyzed vortex structures in a 
polarized Fermi gas at unitarity. The lowest bound state 
will be visible via phase-contrast imaging, together with 
a quantum phase transition at a critical spin polarization. 
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